Fractional photon-assisted tunneling for Bose-Einstein condensates in a double well 
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Half-integer photon-resonances in a periodically shaken double well are investigated on the level 
of the A'^-particle quantum dynamics. Contrary to non-linear mean-field equations, the linear N- 
particle Schrodinger equation does not contain any non-linearity which could be the origin of such 
resonances. Nevertheless, analytic calculations on the A'^-particle level explain why such resonances 
can be observed even for particle numbers as low as A'^ = 2. These calculations also demonstrate 
why fractional photon resonances are not restricted to half-integer values. 
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I. INTRODUCTION 

Tunneling control of ultra-cold atoms via time-periodic 
shaking [11 12 [31 H] of potentials is currently established 
as an experimental method both on the single particle 
level |S| and on the level of Bose-Einstein condensates 
(BECs) An interesting effect is an analog of photon- 
assisted tunneling in periodically shaken systems of ultra- 
cold atoms. It was predicted theoretically both for the 
case that the driving frequency matches the potential dif- 
ference between neighboring wells [H |7] and for the case 
that the driving frequency is resonant with the interac- 
tion energy 8 . The rt-photon resonances essentially are 
a single-particle effect which survives interactions; one- 
and two-photon resonances have been observed experi- 
mentally for BECs in periodically shaken lattices . The 
"photons" are time-periodic potential modulations in the 
kilo-Hertz-regime . 

However, photon-assisted tunneling is not restricted 
to integer photon resonances. Also half-integer Shapiro- 
like |10| resonances have been predicted numerically both 
on the mean- field ( Gross-Pit aevskii) level and on the level 
of the multi-particle quantum dynamics (down to A'^ = 2 
particles) |3j- While the occurrence of higher or lower 
harmonics in non-linear equations is easy to understand 
qualitatively, it is not clear a priori how these resonances 
should occur in the linear A^-particle Schrodinger equa- 
tion. Thus, analytic calculations which can explain the 
occurrence of such resonances within the linear quan- 
tum dynamics will explain how effective non-linearities 
can arise from linear dynamics even for small particle 
numbers. Realistic experimental values for the number 
of atoms in a double well can be of the order of 1000 
atoms [11] for BECs and down to less than 6 atoms [12] 
for few-atoms experiments. 

Often, Floquet-states [T^] help to understand the 
physics of BECs in periodically driven systems [HI [151 
[T51 [Tf] . The focus of the present paper, hes on a dif- 
ferent approach: analytic calculations on the A^-particle 
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level developed in Ref. [H] (cf. Ref. [Ej). By assuming 
the experimentally realistic initial condition of all par- 
ticles being in one well [TT, the calculations are done 
analogously to the time-dependent perturbation theory. 

The paper is organized as follows: after introducing the 
two- mode model for a BEC in a double well (Sec.|ll]), we 
develop the technique to calculate half-integer resonances 
in Sec. |III[ A crucial test is to show that the analytic 
result vanishes in the limit of non-interacting particles 

(Sec, [rvt 

Sec. [VI 



Other fractional resonances are discussed in 



II. THE MODEL: A BEC IN A DOUBLE WELL 

Bose-Einstein condensates in double-well potentials 
are interesting both experimentally and theoretically [ll] 
[20l[2ll[22l[23l[ai[25l[26]. In order to describe a BEC 
in a double well, we use a model originally developed 
in nuclear physics j27j : a multi-particle Hamiltonian in 
two- mode approximation [28j . 



H 



/i(^o + Misin(tJt)) (^4^2 -clci) , (1) 



where the operator annihilates (creates) a boson in 
well j; Ml is the tunneling splitting, K^lq is the tilt be- 
tween well 1 and well 2 and /i/ii is the driving amplitude. 
The interaction between a pair of particles in the same 
well is denoted by 2/ik. 

The Gross-Pitaevskii dynamics can be mapped to that 
of a nonrigid pendulum |20j . Including the term describ- 
ing the periodic shaking, the Hamilton function is given 
by: 



Nk 



— \J\ — cos(0) 



(2) 



where and z are canonically conjugate variables. The 
quantity z/2 is the population imbalance with z/2 = 0.5 
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[zjl = —0.5) referring to the situation with all particles 
in well 1 (well 2). The corresponding observable on the 
A''-particle level is given by 



N 



2N 



(3) 



For integer photon-assisted tunneling, the potential 
difference between both wells, 2h^o has to be bridged 
by an integer number of photons: 

2hfio — nhoj , 71 = 1, 2, . . . (4) 

The 1/2-integer resonance occurs for 

2^A*o = ]^fv^ ; (5) 

for an interacting bose gas, these resonances are further- 
more shifted [21 • 

For some parameter regimes (especially for interac- 
tions comparable to the onset of the self-trapping tran- 
sition [H]), the differences between mean- field (Gross- 
Pitaevskii) dynamics and the iV-particle quantum dy- 
namics can be quite remarkable ,2-9] • However, when 
concentrating on the (experimentally measurable 
time-averaged population imbalance. 



{J: 



N 



1 

f 



N 



dt. 



(6) 



for low interactions, the qualitative agreement between 
mean-field and A^-particle dynamics for the occurrence of 
both integer and half-integer photon-assisted tunneling is 
excellent [3 . 

Photon-assisted tunneling is clearly visible in the ex- 
perimentally measurable time-averaged population im- 
balance (|6|. Figure [l] shows integer resonances (|4]), 
namely the one-photon peak with uj « Sfi and the two- 
photon peak at a; ~ 1.5f2. Furthermore, there are pro- 
nounced fractional-integer resonances at w w 651, w 2ri 
and oj « 1.2ri corresponding to the 1/2, 3/2 and 5/2- 
photon peaks. While some of the resonances disappear [31 
for specific choices of the driving amplitude, the initial 
phase of the periodic driving (cf . ) does not influence 
the occurrence of resonances in the situation investigated 
in this manuscript. 

Figure |2] shows that it is not essential to start with all 
particles in one well in order to observe photon-assisted 
tunneling. Both for the ground-state of the untilted, un- 
driven system (for which the initial population imbalance 
is zero) and for the ground-state of the tilted system with 
an initial population imbalance of ~ 0.467, the main res- 
onances of Fig. [T] where all particles were initially in 
well 1, can easily be identified. 

Figure [3] displays the half-integer resonance for N = 2 
particles. Contrary to what was observed for both larger 
particle numbers and for mean-field, the position of the 
1/2 photon resonance does not shift with increasing en- 
ergy. A first test of our analytic calculations towards 
the end of the next section will thus be to explain this 
feature. 
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FIG. 1: (Color online) Time-averaged population imbalance 
{Jz)t/N with averaging time T = 100/f2, static tilt 2fio/Q = 
3 and interaction Nk/Q = 0.4 for N = 2 particles initially in 
well 1. The red (full) line with driving amplitude 2fj,\/u = 1.8 
(close to the maximum of the Ji-Besselfunction) has, among 
others, a pronounced resonance at the driving frequency lo = 
30 (one-photon resonance) and also a weaker resonance at 
UJ — 6n (l/2-photon resonance). As the black (dashed) line 
shows, these resonances are also visible when the system is 
driven with an initial 7r/2 phase shift (^i cos(tjt) instead of 
fii sin(ci;t) in Eq. ([ij), while they are strongly suppressed with 
driving amplitude 2^\/uj = 3.83 (close to the first zero of Ji) 
illustrated by the blue (dotted) line. 




FIG. 2: (Color online) Time-averaged population imbalance 
{Jz)t/N with averaging time T = lOO/fi, static tilt 2^o/Q, = 
3 and interaction Nn/Q, — 0.4 with N — 2 particles for two 
different initial states. The dotted line corresponds to the 
ground state of the tilted system, while the full line displays 
{Jz)t/N with the ground state of the untilted system as initial 
state. In both cases the 1/2-photon resonance at ui = 6Q 
appears. 



III. ANALYTIC CALCULATIONS 

In order to analytically describe the time-evolution of 
the interacting system, the Fock basis {v) = |iV — ly,!^) 
is used. The label — ... TV refers to a state with 
N — V particles in well 1, and v particles in well 2. The 
Hamiltonian ([T]) now is the sum of two (7V-t- 1) x {N +\)- 
matrices, 



H = Hoit)+H, 



(7) 
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the A^-particle Schrodinger equation governed by the 
Hamiltonian ([T|): 

ihd^{t) = {iy\Hi\iy+l)K{t)a^+i{t) 

+ {iy\H,\j^-l)K_i{tya,^i{t) . (10) 



In Eq. (lOl, the notation a^i{t) = aN+i{t) = 0, was 
used; the phase factors are given by: 

h^,{t) = exp {i [2{N -l-2v)Kt + 2/ioi - 2/ii cos(wi) /uj]) 

(11) 

with — cos(ti't)/a; = ^^sm{ijjt')dt' . To simphfy the ex- 
pression for subsequent integrals, one can use the expan- 
sion in terms of Bessel functions 30 




FIG. 3: (Color online) Time-averaged population imbalance 
for N = 2 particles in a periodically shaken double-well po- 
tential as a function of both driving frequency ui and driving 
amplitude /ii for a static tilt of 2jj,o/Q — 3. Upper panel: 
Nk/Q, = 0.2; lower panel: Nn/Q — 0.4. The averaging time 
is QT = 100; the values plotted are shifted such that corre- 
sponds to all particles having always stayed in the first well (a 
value of 1 would correspond to all particles being in well 2). 
Surprisingly, contrary to the case of A'' > 2 or the mean-field 
case [3] , for N — 2 particles the resonance does not shift with 
increasing interactions. 



While the non-diagonal matrix Hi is given by the 
tunneling-terms of Eq. ([T|), the diagonal matrix Hq in- 
cludes both the interaction between the particles and 
the applied potential difference. For the solution of the 
Schrodinger equation 



the ansatz 



{v\ip{t)) = a^{t)exp 



{,y\Hoit')\iy)dt 



(8) 



(9) 



turned out to be useful [18] . 

Within this framework, a set of differential equations 
was derived [18 which is mathematically equivalent to 



s{cijt) 



- J2 J;c(^)*'e^'^ 



k — — OCi 



Equation (10 1 furthermore needs 



{,y\Hi\n)= - 



2 

nn 



Si^m+WN - n\/n + 1 
5u,n-WN - n + \^/n , 



(12) 



(13) 



where Sn^m is the Kronecker delta (which is zero ex- 
cept for n = m where Sn,n = 1)- The idea is to pro- 
ceed along the lines of time-dependent perturbation the- 
ory |31j . Starting with a typical experimental initial con- 
dition such that all particles are in the first well [11], one 
has in zeroth order perturbation theory: 



1. 



0, 



(14) 



where = J2k=o 
one gets: 



(fc) 



In first order perturbation theory, 
aW(t) = (15) 



ii v 1 and 



(16) 



Using Eqs. ([rT|) and ([121 one thus has 
,(1) 



al'-'it) ^i'-^VN V i^Jk{2^ii/uj) i exp{i(Tkt')dt' 



with 



(jfc = few - 2/10 - 2(A^ - 1)k 



(17) 



(18) 



Therefore, after solving the integral, Eq. ( [17| ) is a sum of 
time-periodic functions except for the special case with 
CTfe = which recovers the integer photon resonances of 
Eq. ^ investigated in Refs. [3j [9]. While for a double 
well as in Ref. [3 , the population imbalance is ideal to 
investigate photon-assisted tunneling, the experiment [S] 
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was performed in an optical lattice. The signatures of 
photon-assisted tunneling were seen in the width of the 
BEC after expansion in the shaken lattice. Surprisingly, a 
dependence was measured. While this might be inter- 
preted as being an indication for transition from ballistic 
to diffusive transport |9] , the present experiments cannot 
exclude other explanations. The J^-dependence could ei- 
ther be an interaction-induced effect [32 or the result of 
an effective average over the precise instant within the 
cycle at which the current is measured [4]. 

As the aim of the present paper is to understand the 
fractional photon peaks like the interaction induced half- 
integer resonances of Ref . ^ , we can discard the integer- 
photon resonances characterized via ak = [33] and thus 
write 



0.01 



0.005 -; 




a['\t)^t^^ V ^'^Jf,{2^,,/uJ)^—^ (19) 

fc— — oo 



In second order perturbation theory one has (see the ap- 
pendix lAl): 



,(2) 



FIG. 4: (Color online) Time-averaged probability (averaged 
over TO = 10) that two particles have tunneled to the other 
well for the 1/2-photon resonance (A'^ = 2, fio = 30/2, 
u — 60, K — 0.2O). Wide (green/grey) line: num erica l data; 
dashed (black) line: perturbation theory (cf. Eq. (A8l). The 
probabilities displayed here should be measurable experimen- 
tally (see Ref. [12]). 



(t) = ^i^^ Vn^/N - l^/2 (20) where the expression 

oo oo 

X ^ ^ i'=iVfe(2^i/w)j£(2^i/w) 

A;— — OO ^— — oo 

exp{i(Tkt') - 1 



(^2^^ (i))oscil 



contains oscillatory 
terms which can be found in Eq. (A7). The convergence 



of this sum is ensured both by the scaling of (Jk on k and 
the behavior of Bessel functions with increasing k [50] 



/ 



Jk{z) 




(25) 



with 



ai = luj- 2/io - 2(iV - 3)k 



(21) 



Again, ai — Q can be discarded because it corresponds to 
integer photon resonances. However, if -I- cr^; = then 

(2) 

02 does have parts which increase linearly in time. In 
order to see that this indeed corresponds to a half-integer 
resonance, we choose N — 2 and llj/2 = 2^q. This implies 
CTfc = kij — 10 12 — 2k and 5^ = luj — uj 12 ^ 2k and the 
condition 



CTf -I- (Tfc 







(22) 



thus becomes independent of the interaction; it results in 
the simple equation 



(23) 



The above reasoning explains why we observe no shift 
of the resonance with increasing interaction in the nu- 
merics displayed in Fig. [3] The amplitude to find both 
particles in well two is given by (appendix [A|: 



combined with the fact that J^k{x) = {—!)'' Jk{x) (for 
integer k) . Figure [4] shows good qualitative agreement 
between the analytic and numeric calculations for the 
time-averaged probability that both particles, which ini- 
tially have been in the first well, have tunneled to the 
second well. Already perturbation theory in the first or- 
der, in which the half-integer resonance becomes visible, 
correctly describes the occurrence of maxima and min- 
ima in the probability for both particles to occupy the 
second well. 



IV. HALF-INTEGER RESONANCES 
DISAPPEAR IN THE LIMIT OF LOW 
INTERACTIONS 

Despite the agreement displayed in Fig. |4| at the first 
glance Eq. (24 1 seems to contain a flaw: numerically. 



4"w = (4'^w) 



oscil 



(24) 



f^2 °o ^ 

-^fc (2/ii/w) Ji-fc (2^i/w) — 



we observe that the half-integer resonance disappears for 
zero interaction. However, there seems to be a sum of 
non-zero terms proportional to t even for k = 0. As it 
is not obvious that these terms cancel, the next step will 

(2) 

be to demonstrate that 03 indeed approaches zero for 
vanishing interaction. 

As shown in the appendix, the terms pr opo rtional to t 



(2) 

in 03 are due to situations such that Eq. ( 22 1 is fulfilled 



1 



k — ~OD 



In the limit k — s- this results again in the condition ( 23 1 
independent of the particle number. The part of 
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which increases hnearly in time is thus proportional to 

oo ^ 

A= y] Jfc(2/^i/w)Ji_fe(2^i/a;)— . (26) 

• * /Ti. 



Dividing the sum into two parts (X^fcLi ' ' ' + J2e=-oo • ■ ■ ) 
and then setting 1 — £ — k, one obtains 
oo . ^ 

A = Y, Jfe(2Aii/c^) Ji-fc(2/ii/cc;) - 



k=l 



(27) 



In the limit k ^ 0, the position of the half-integer reso- 
nance approaches the value for N = 2 particles. There- 
fore, one has ak — kcu — lu/2 and thus ai-k = ^(^k which 
implies 



A = Q 



(28) 



Thus, in agreement with the numerics, the half-integer 
photon peak disappears with vanishing interactions. 



V. FRACTIONAL INTEGER RESONANCES 

Fractional integer resonances are not, however, re- 
stricted to the half-integer resonances investigated nu- 
merically in Ref. [31 and analytically in Sec. IV For 
= 3 particles and a driving frequency such that 
cij/3 = 2/io, the condition 







(29) 



with 



CT„ = mw-2^o-2(7V-5)K (30) 
(throughout this section: iV = 3) is fulfilled for 



k + £ + m=l 



(31) 



The amplitude to find three particles in well 2 again 
contains oscillatory terms, the term which becomes the 
leading-order term for large t can be obtained by a cal- 
culation analogously to the half-integer resonance in ap- 
pendix |A] 

3173 



(4") 



linear 



Jfc(2Mi/^) 



(32) 



X J^(2^i/w)Ji_£_fc(2Mi/w) 



0'fe(CTfe + CTf) 

This one-third photon resonance can indeed be observed 
in the numerics (see Fig. [s]). As this resonance only 
occurs in third order perturbation theory (rather than 
second order for the half-integer resonances), the ampli- 
tudes would be rather small for interactions as in Fig. |3] 
However, choosing an also realistic value of Nn/il = 1.5 
leads to a time-averaged population imbalance with a 
peak-height of the same order of magnitude as in Fig |3] 
In a similar manner, smaller fractions could be treated 
in higher order perturbation theory. As the resonances 
thus are a higher order effect, they will tend to decrease. 




0.15 



2 2\i^/()) 



FIG. 5: The plot illustrates the 1/3-resonance for TV = 3 
particles and interaction parameter Nn/Q — 1.5. Initially, all 
three particles were in the lower well; the static tilt is again 
given by 2po/^^ = 3. The time averaged probability (averaged 
over time Til = 100) to find all particles in the upper well 
as a function of the driving frequency uj/Q, and the driving 
amplitude 2^i/uj has a clear peak at w ~ 



VI. CONCLUSION 



Contrary to the integer-photon peaks [5], fractional- 
integer photon peaks cannot be explained by simply 
replacing the time-dependent Hamiltonian by a time- 
independent Hamiltonian with renormalized tunneling 
frequencies. As half-integer resonances already appear 
for two particles in a double well, this experimentally rel- 
evant case [12] was investigated both numerically and an- 
alytically. The perturbation calculations can explain for 
which parameters the non-integer resonances occur. As 
the fractional-integer resonances are only visible for finite 
interactions between the particles, they allow to investi- 
gate beyond single-particle effects for very small particle 
numbers. Experiments similar to Ref. [12j could thus ver- 
ify fractional-integer peaks in photon assisted tunneling 
and thus help to understand the emergence of effects sim- 
ilar to the non-linearities of a mean-field approach well 
below the limit N ^ oo. 
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APPENDIX A: SECOND ORDER 
PERTURBATION THEORY 



with 



When solving the integral 



Ik.e = I - — [exp (i((Tfc + ae)t') - exp {iaet')] dt' (Al) 



in Eq. ( 20 1 , one can again assume CTfe and CT| as 



CTfc = and ai — would correspond to the inter-photon 
resonances discarded here. It then remains to distinguish 
cases with 



<Jk + <7e = , 



(A2) 



which turn out to be the origin of the half-integer reso- 
nance, from those for which this equation is not fulfilled. 



If Eq. (A2) is fulfilled, one has 



i(jk 



t H (exp(— i(T/jt) — 1) 



(A3) 



otherwise 



1 



exp[z(tTfe -|- ai)t] — 1 exY){iuit) — 1 



(A4) 

Collecting all terms given by Eq. ( |A3[ ), one has the 
leading-order contribution: 



leading— order 



(A5) 



Jk{x) Ji-k{x) 



k — -~c>o 



CTfc 



exp(— icTfet) — 1 



lak 



X = 2/ii/(jj 



(A6) 



and an oscillatory part 



(af(t)) 



(A7) 



E E 



Jk{x)Ji{x)i 



k+l-l 



k=-oo e^l-k 

exp(z((Tfc + ffe)t) — 1 exp(iaet) — 1 



{i{ak + (Ti)) 



While Eq. (A5l includes the leading-order behavior 



for large times and most parameters, it vanishes in the 
limit /ii — !■ 0. Thus to evaluate the analytic formula with 
the help of a computer algebra program, we include the 
only non-vanishing term for /xi = to obtain the data 
displayed in Fig. [4] 



\ / approx \ / 



leading— order 

Joix)Joix)i-^ 



(A8) 



exp(i(o-o -I- q-q)^) 
(i(cro + S^o)) 



exp^icTot) — 1 
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